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Long-lifetime Polariton BEC in a Microcavity with an Embedded Quantum Well and 

Graphene 
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We study the propagation of a Bose-Einstein condensate (BEC) of long-lifetime exciton polaritons 
in a high-quality microcavity with an embedded semiconductor quantum well or a graphene layer 
using the Gross-Pitaevskii equation. It is shown that in these heterostructures the BEC of the long- 
lifetime polaritons can propagate over the distance up to 0.5 mm. The obtained results are consistent 
with the recent experimental observations for GaAs/AlGaAs microcavity. It is demonstrated that 
the BEC density in a polariton trace in a microcavity with embedded graphene at large distances 
from the excitation spot is higher for the microcavity with higher dielectric constant. It is also 
predicted that the propagation of a polariton BEC in a microcavity with graphene is dynamically 
tunable by changing the gap energy, that makes it potentially useful for applications in integrated 
optical circuits. 

PACS numbers: 71.36.-l-c, 71.35.Lk, 73.21. Fg, 78.67.Wj 
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I. INTRODUCTION 

Polaritons are the quantum superposition ofphotons 
and excitons in semiconductors or graphene [l|-l3|. In 
the past years, formation and spreading of the po- 
lariton Bose-Einstein condensate (BEC) in quasi-two- 
dimensional structures attract attention of experimen- 
talists and theoreticians from the point of view of funda- 
mental physics at nanoscales and due to potential ap- 
phcations of polaritons in integrated circuits for opti- 
cal and quantum computing [j-[6|. Because the effective 
mass of the polaritons is much smaller than the atomic 
masses, the BEC transition temperature is much higher 
than that for atomic BECs [7|. Another important fea- 
tures of the polariton BEC dynamics is that it is essen- 
tially non-equilibrium due to the finite polariton lifetime. 
The polariton decay is mainly caused by the leakage of 
photons from an optical microcavity and usually the po- 
lariton lifetime is of the order of a few picoseconds, de- 
pending on the quality Q-factor of the microcavity [3|. 
It was established that the polariton BEC can move in 
a planar microcavity in response of the external force or 
field [sl, Q. However, in "traditional" experiments the 
characteristic length over which the BEC can spread in 
the microcavity is limited due to relatively short lifetime 
of the polaritons [9| . 

The signatures of the polariton BEC were first ob- 
served in GaAs-based microcavity |1Q|-[15| and then in 
CdTe-based [l6| and GaN-based microcavities |17j . 
There are two advantages for observation of a polariton 
BEC in such heterostructures over a conventional exciton 
BEC [l8|. First, polaritons have an effective mass, which 
is ^ 10~^ times lower than the exciton mass, that results 
in much higher critical temperature of the BEC transi- 
tion for polaritons than that for the exciton BEC at the 
same particle density [l|, [l9|, [20| . Second, the polariton 
BEC is less sensitive, compared to the exciton BEC, to 
the disorder caused by unavoidable crystal defects and 
impurities. It is worth noting that the presence of the 
defects results in the exciton localization in a fluctuating 



potential in a sample and can destroy an exciton BEC 

It was found [IMS, E E [UHlij that for optimal 
conditions of the polariton BEC formation the quality Q- 
factor of the optical microcavity should be enough high. 
The increase of the microcavity Q-factor leads to the in- 
crease of the polariton densities in the microcavity and, 
also, to longer polariton lifetime thus presenting more 
favorable conditions for the BEC observations. Also, 
enough strong exciton-photon coupling (Rabi splitting) 
is needed thus, requiring large quantum well (QW) exci- 
ton oscillator strength. Enough large polariton-polariton 
interaction scattering cross-section is also an important 
factor, which provides efficient thermalization in the po- 
lariton system and hence short enough relaxation time 
for the transition of the system to a BEC (see Ref. [19| 
for the discussion). Rapid expansion of a spatial coher- 
ence and the polariton BEC formation dynamics in a 
GaAs-based microcavity has been recently observed in 
Ref. [30J. 

In the present Article we consider the dynamics of 
polaritons formed by the cavity photons and the exci- 
tons in GaAs a quantum well and graphene embedded 
in a microcavity. Our studies are motivated by experi- 
ments [8] and [31] where the directional propagation of 
exciton polaritons have been observed in the presence of 
an external force. In particular, in Ref. [31| the possi- 
bility of synthesis of a semiconductor microcavity with 
the Q-factor exceeding 10^ has been demonstrated for 
the first time. That allowed one to observe the po- 
laritons with extremely long lifetime ~ 100 ps, that is 
^ 30x longer than that in the traditional experiments 
[35. It was found [3i| that the long-lifetime polaritons in 
a planar microcavity could propagate over a macroscopic 
scale of a fraction of millimeter order. It was supposed 
in Ref. [31| that a BEC coherent flow formed at high 
enough densities of the long-lifetime polaritons have pro- 
moted their long-distance propagation. In this Article, 
we demonstrate that the polariton dynamics observed in 



Ref. [3l| is consistent with that predicted from the Gross- 
Pitaevskii equation for a polariton BEG. Specificahy, we 
show that under the action of an external force, a long- 
hfetime polariton BEG can propagate over the distance 
~ 0.5 mm from the excitation spot. It is demonstrated 
that the polariton-exciton interactions can also signifi- 
cantly change the spatial distribution of the polaritons 
in the BEG. Based on the obtained results, we propose 
an observation of long-lifetime polariton BEG formed in 
a graphene layer embedded into a high-Q microcavity. It 
is shown that for graphene there is an additional control- 
ling parameter - the energy gap in the electron and hole 
excitation spectra - that allows one to govern the polari- 
ton BEG propagation. The possibility to dynamically 
adjust the properties of the polariton BEG by changing 
the gap, for example through application of the electric 
field, makes a microcavity with embedded graphene a 
promising candidate for the use in adaptive optical cir- 
cuits. 

Our article is organized in the following way. In Sec. 
II the dynamics of a Bose-Einstein condensate of long- 
lifetime exciton polaritons in a microcavity with an em- 
bedded semiconductor quantum well or a grapheme layer 
is considered based on the Gross-Pitaevskii equation. 
The details of the numerical simulations are presented 
in Sec. III. Sec. IV presents the results of the simulations 
of the polarition BEG propagation in a microcavity with 
an embedded quantum well and an embedded graphene 
layer. Here, the effect of the polariton lifetime on the 
BEG dynamics is discussed and the comparison of the 
BEG propagation in a microcavity with an embedded 
quantum well and graphene is presented. Finally, the 
conclusions follow in Sec V. 

II. DYNAMICS OF AN EXCITON-POLARITON 
BEG IN A MICROCAVITY 

At temperature much lower than the BEG transition 
temperature the dynamics of the exciton polariton BEG 
is described by the Gross-Pitaevskii equation 
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+^^(r,t)|^(r,t)p-zn7^(r,t) + inG(r). (1) 

In Eq. ([!]) ^(r,t) is the condensate wave function which 
depends on the coordinates r = (x^y) in the plane of 
the microcavity and time t, m is the polariton mass, g 
is the effective interaction strength proportional to the 
polariton-polariton scattering amplitude, 7 = l/2r is the 
decay rate, r is the polariton lifetime, and G{r) is the 
source term. The term U{r) = /7ex-poi(^) + ^ext(^) de- 
scribes the repulsive interaction /7ex-poi(^) of the polari- 
tons with those excitons in a quantum well, which are not 
coupled with the microcavity photons, and the potential 
energy due to the external force, Uext{r). Importance of 
the interactions of polaritons with the exciton cloud has 
been recently discussed in Ref. [32|- The corresponding 
interaction energy is 



where Xq = Xk=o is the Hopfield coefficient for the ex- 
citon component in the polariton wave function taken 
at the wave vector /c = 0, ^ex-ex is the exciton-exciton 
repulsive interaction constant and nex(^) is the exciton 
density in the cloud. The Hopfield coefficient for arbi- 
trary wave- vector k is 
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where D,r is the Rabi sphtting, 
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is the energy of the lower polariton branch, Ex{k) is the 
exciton energy and Ec{k) is the cavity photon energy as 
functions of k [33]. Because the polaritons in the BEG 
are condensed into the state with small wave- vectors /c, 
the Hopfield coefficient in Eq. (J2]) is taken at /c = 0. We 
approximate the exciton density in the cloud via a Gaus- 
sian profile as nex(^) = noexp(— It — rop/a^) where no is 
the density at the center of the cloud and a is the charac- 
teristic size of the cloud. The Gaussian profile gives the 
required bell-like dependence for the exciton cloud den- 
sity in the excitation area |r — ro| ^ a and also captures 
the sharp drop of the exciton density at large distances 
|r — ro| ^ a [7]. Thus, the polariton-exciton interaction 
(J2]) is represented as 



^ex-poi(^) = Uoexp{-\r -rof/a'^), 



(5) 



where Uo = Xg^ex-ex^o- To study the effects of the 
polariton-exciton cloud interaction, we vary the interac- 
tion strength Uq as described below. 

It has been demonstrated in recent experiments [8|, [31| 
that the microcavity polaritons can be accelerated by an 
external force due to the spatial gradient of thickness 
of the microcavity. In this approach, the average force 
acting upon a polariton wave packet at the point r is 
F{r) = —\/Ec{r)^ where Ecir) is energy of the polari- 
ton band taken at the in-plane vector of t he p olariton 
/c = [8]. In particular, in experiments [sl, [3l| the mi- 
crocavity width, and hence the energy Ec{r)^ has been 
a linear function a spatial coordinate (the wedge-shape 
geometry of the microcavity) that resulted in a constant 
force F{r) = Fq acting upon the polaritons. To describe 
the effect of the force, we set the corresponding potential 
energy as 



Ue^t{r) = -Fox. 



(6) 



^ex-pol(^) = X^g, 



ex— ex'f'ex 



(r), 
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To explore the polariton BEG dynamics in a macro- 
scopically large sample we consider the case where the 
source of polaritons, G{r), is localized in space [29|]. 
Within the used approximation, the source term is rep- 
resented as G{r) = Go exp(— |r — rop/a^) with the same 
characteristic size a as for the exciton cloud size. 
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The polariton-polariton interaction strength is [33|, |3^ 
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FIG. 1: Time dependence of the total number of polaritons in 
the BEG for Fq = 13 meV/mm and Uo = 3 meV for different 
pumping rates Co. 

Fohowing the conditions of the experiments [sl, l31j , 
we varied Uo from to 10 meV and Fq from to 
13 meV/mm. For the long-hfetime polaritons in a micro- 
cavity with the quality factor Q ~ 10^ with an embedded 
semiconductor QW we set r = 100 ps |31], m = 10~^mo 
where mo is the free electron mass [3], and g = 10~^ 
meV/im^ [9]. 

Formation of the polariton BEC in 
has recently been predicted in Ref. [3; 
effective mass in graphene is 



apped graphene 
. The polariton 



m = 2(m^^ + cLcl^feixK) 



(7) 



where ?7iex is the exciton effective mass in graphene, L^ is 
the length of the microcavity, e is the dielectric constant 
of the microcavity, and c is the speed of light in vacuum. 
In the present work we consider the case of zero detuning 
where the cavity photons and the excitons in graphene 
are in resonance at /c = 0. In this case the length of the 
microcavity is chosen as |35| 



Lr = 
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{2S-Vo^C/S^)^' 



(8) 



where Vb = e'^/Aneeor' , C = {Hcvf)'^ /Sneeor'^^ S is the 
gap energy in the quasiparticle spectrum, and vf is the 
Fermi- velocity of the electrons in graphene. The param- 
eter r' is found from the following equation. 



26^{2S-hujy^ 
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where D = e'^/ATreeo. For dipolar excitons in 
GaAs/AlGaAs coupled quantum wells, the energy of the 
recombination peak is huj = 1.61 eV [36|. We expect simi- 
lar photon energies in graphene. However, its exact value 
depends on the graphene dielectric environment and sub- 
strate properties though. The exciton effective mass in 
Eq. ([7]) is given by 
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where a^ = 2i\eosT? jvnre^ is the two-dimensional Bohr 
radius of the exciton, vfir is the exciton reduced mass. 

Formation of the excitons in graphene requires a gap 
in the electron and hole excitation spectra, which can be 
created, for example, by external electric field |37H39|. 
We do not restrict our consideration to a specific mecha- 
nism of the gap formation. As follows from Eqs. ([7|)- (pT]) , 
both the polariton mass m and the interaction strength 
g in gapped graphene depend on the gap energy ^. In 
our studies, we varied ^ from 0.1 eV to 0.5 eV that are 
the representative values for graphene. The polariton 
lifetime is mainly defined by the leakage of the photons 
through the mirrors and only weakly depends on the exci- 
ton lifetime [3|,[31|. Therefore, we use the same polariton 
lifetime r = 100 ps for graphene and for a semiconductor 
quantum well. 



III. NUMERICAL SIMULATIONS 



To study the polariton BEG dynamics, we numeri- 
cally integrated Eq. ([1]). At the initial moment t = 
we set ^(r,0) = 0. As is demonstrated in Sec. IIVI be- 
low, the long-lifetime polaritons accelerated by the exter- 
nal force can propagate over distances up to /^ax '^0.5 
mm from the excitation spot. The characteristic wave- 
length of the polariton wave packet at the distance 
/ ^ / max can be estima ted quasiclassically as Amin ^ 
27rfi/Y^2m(Fo/max + ^o) that gives Amin ~ 9.5 x 10~^m 
for Fo = 13 meV/mm and [/q = 10 meV. At shorter 
distances / < /max one has A > Amin- To resolve the 
oscillations of the polariton wavefunction with the wave- 
lengths A > Amin, we choose the numerical grid spacing 
Ax = 1.5 X 10~^ m. The dimensions of the numerical 
box are L^ x Ly = 6000 x 3000 (Ax)^ The long side of 
the numerical box is arranged parallel to the direction of 
the constant force Fq. We use the periodical boundary 
conditions at the numerical box boundaries. The unit of 
time is equal to At = 2m{/^xf /h = 3.9 x lO"^"^ s. We 
numerically integrated Eq. ([T]) over time with the 4*^ or- 
der Runge-Kutta scheme at the time step h = 10~^At. 
An explicit 4*^ order accuracy method has been used to 
represent the Laplacian operator in Eq. ([1]). We took 
advantage of the graphical processing unit GUDA pro- 
gramming [40] that allowed us to achieve ~ 20 x increase 
of the performance with respect to the shared-memory 
parallel version of the code. This allowed us to perform 
the simulations for this macroscopic system within prac- 
tically realizable time intervals. The overall computation 
time for a single run was ~ 10 hours on NVIDIA Tesla 
S2050 card. 



IV. PROPAGATION OF A BEG OF 

LONG-LIFETIME POLARTIONS IN A 

MIGROGAVITY 

A. Polartion BEG in a microcavity with an 
embedded quantum w^ell 

In this Section we focus on the propagation of a BEG 
of long-lifetime polaritons in a microcavity with an em- 
bedded semiconductor QW. Fig. [1] shows the time de- 
pendence of the total number of polaritons Np{t) = 
/(i^r|^(r,t)p after the source it turned on. It follows 
from Fig. [1] that after the system approaches a steady 
state at t ~ 6 X 10^ At ^ 230 ps, for the pumping rate 
Go = 0.1Ax~^At~^ the total number of polaritons is 
Np ^ 5.3 X 10^ that is a representative value for exper- 
iments with exciton polaritons in a semiconductor mi- 
crocavity [311 . Below we used this pumping rate in our 
simulations. 

Propagation of long-lifetime polariton BEG in the 
quantum well plane is shown in Fig. [2l The center of 
the excitation spot is positioned at ro = (xo,^o) where 
xo = yo = 1500Ax. The evolution of the long-lifetime 
polariton BEG distribution with increasing of the force 
Fq is given in Fig.[2^-c. Fig. [2^ shows that in the absence 
of the force the polariton spot is symmetric in the (x, y) 
plane and has a characteristic size of ^ 500Ax = 75 /im 
that is, ~ 7.5 a. With the increase of the force a polariton 
"trace" is formed and at Fq = 13 meV/mm the length 
of the trace reaches ^ 3.5 x 10^ Ax ^ 500 /am (Fig. [2fc). 
To characterize the steady-state polariton distribution in 
details we plot in Fig. [2li the polariton BEG density at 
y = yo that is, a cross-section of Fig. [2^-c along the line 
passing through the center of the excitation spot. The in- 
set in Fig. [2]i shows that for the force Fq = 13 meV/mm, 
the maximum of the polariton density is shifted from the 
center of the excitation spot xq to x ^ 1620Ax. 

To prove that the shift of the density maximum is 
caused by the interaction of the polaritons with the ex- 
citon cloud, we investigated evolution of the polariton 
density with changes of the maximum effective interac- 
tion strength Uq. Fig. [3^ demonstrates changes in the 
dimensionless polariton density along the line y = yo 
with the rise of Uq from 1 meV through 3 meV. It is 
seen that the position of the maximum of the density is 
shifted from x ~ 1500Ax at [/q = 1 meV to x ~ 1620Ax 
at [/o = 3 meV. Fig. [Sb-d shows the contour lines for 
the polariton density n{r) for different Uq. From Fig. [Sb 
it is evident that at [/q = 2 meV and y = yo = 1500 Ax 
the spatial polariton distribution has two maxima posi- 
tioned at X ~ 1500Ax and 1620Ax. While increasing Uq 
to 2.5 meV the left maximum disappears and the right 
maximum becomes more pronounced, as seen in Fig. [3t. 
With the further increase of Uq, most polaritons are con- 
centrated around the maximum at x ~ 1620 Ax (Fig. 

To clarify the reasons for the shift of the density 
maximum while Uq is increased, we compare in Fig. |4] 
the distribution of the polariton flux in a steady state 



j{r) = (h/2im) (^*(r)V^(r) - ^(r)V^*(r)) for differ- 
ent forces. For Fq = the flux j{r) is distributed sym- 
metrically around the center of the excitation spot. As 
is seen in Fig. H^ the magnitude of the flux, j, reaches 
its maximum at the distance r^ 120 Ax from the exci- 
tation spot center ro- The presence of the maximum 
of j at certain distance from the point ro can be un- 
derstood if one considers the creation and propagation 
of polariton wave packets in the potential energy profile 
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FIG. 2: (a-c) Evolution of a trace of the long-lifetime po- 
lariton BEC with increasing the force Fq for t — 4800At and 
[/o = 3 meV. The colorbars show the dimensionless BEC den- 
sity n(r)Ax^. (d) Evolution of a steady-state polariton den- 
sity at ^ = ^0 for [/o = 3 meV with increasing Fq shown in 
a logarithmic scale and in a linear scale in the region around 
the excitation spot (inset). 
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FIG. 3: (a) Changes in the polariton density at y = yo with 
the increase of Uo from 1 meV to 3 meV for t = 5.8 x 10"^ At. 
(b-d) : Contour plots for the polariton density for Uq — 2 (b) , 
2.5 (c), and 3 meV (d); other parameters are the same as in 
plot (a). The contour lines are labeled with the values of the 



density n(r)Ax^. 
is gray-shaded. 



In plot (b) the area with n(r)Ax > 0.52 



^ex-poi(^)- The polaritons are accelerated due to the 
force Fex-poi(^) = — V/7ex-poi(^) and hence, their veloc- 
ity V gradually rises with the distance from the center. 
At large distances, the velocity v tends to a constant be- 
cause the force Fex-poi(^) vanishes. On the other hand, 
the polariton density n{r) falls at large distances from 
ro due to spreading of the polaritons over the sample. In 
effect, the flux j = nv reaches the maximum at certain 
distance ^ a from the excitation spot center. 
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FIG. 4: Polariton flux j for [/o = 3 meV, t = 5.6 x 10^ At. 
The length of the arrows is proportional to the magnitude of 
the flux |jf|. The point at ro = (1500Ax, 1500Ax) marks the 
center of the pumping spot. 

In the presence of a force Fq , the polariton flux distri- 
bution acquires a strongly asymmetric form. The result 
of the simulations for Fq = 3 meV/mm and /7o = 3 meV 
is shown in Fig. |4]3. The simulations for other forces Fq 
and energies Uq resulted in the polariton flux distribu- 
tions, which were similar to that shown in Fig. |4]3. In 
this case, most of the polaritons move in the direction of 
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FIG. 5: (a) Variation of the density of a long- lifetime polariton 
BEC in a microcavity with embedded graphene with rising the 
gap energy 8. (b) Points: dependence of the polariton BEC 
density at the point r — (2000Ax, 1500Ax) (marked by a 
vertical arrow in plot (a)) with changing 5. Vertical errorbars 
show uncertainty due to small time-dependent oscillations of 
the polariton density about the mean values. The dashed 
curve is a guide to the eye. Fo = 13 meV/mm and [/o = 5 
meV for both plots (a) and (b). 

the force F^. This results in the shift of the maximum 
of the polariton density from the spot center, in agree- 
ment with Fig. [T]i and Fig. [3l It is worth noting that 
at distances r^ a from the excitation spot center there is 
a finite flux directed opposite to the external force Fq. 
However, these polaritons eventually flow around the ex- 
citation spot and, at large distances from the spot, move 
in the direction of the force in agreement with the obser- 
vations in Ref. |31i] . 

B. Polariton BEC dynamics in a microcavity with 
embedded graphene 

Within the same approach, we also studied the propa- 
gation of polaritons in a semiconductor microcavity with 
gapped graphene embedded into it. First, we consider 
the case of the GaAs-based microcavity, for which we 
set e = 13. The effective mass vn of the polaritons and 
the polariton-polariton interaction strength g have been 
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FIG. 6: Dependence of the polariton BEG density in a 
microcavity with embedded graphene at the point r\ — 
(2000Ax, 1500Ax) on the microcavity dielectric constant e. 
The dashed hne is a guide to the eye. 

found from Eqs. ([7|)- (fTT]) for different values of the gap 
energy ^. The results of the simulations for graphene in 
Fig. [5^ shows that the polariton density in the trace grad- 
ually increases with rising the gap energy. To character- 
ize the increase of the polariton BEC density in the trace 
with changing the gap energy we present in Fig. [SJd the 
dependence of the dimensionless polariton density at the 
point r\ = (2000Ax, 1500Ax) on the gap energy 6. The 
point ri is positioned at the distance of 500Ax ~ 75 /im 
from the excitation spot center vq. The position of the 
point ri is marked by a vertical arrow in Fig. [5^. Fig. [SJd 
shows that the polariton density at the point ri increases 
in 10. 2x, from n ^ O.OISA^-^ to n ^ 0.133Ax-^ with 
changing 5 from 0.1 eV to 0.5 eV. Additionally, the shift 
of the density maximum in the direction of the external 
force Fq is more pronounced for larger S^ as shown in 
Fig. [5^. 

The microcavity can be synthesized from different 
semiconductor materials. To study the effect of the mi- 
crocavity material we investigated the propagation of the 
long-lifetime polariton BEC in graphene in microcavities 
formed by semiconductor materials CuBr, ZnSe, CdTe, 
and ZnO [41| that is, by those materials, which were 
utilized for the experimental studies of microcavity po- 
laritons [11, [13, El, |42|-[i3 • According to Eqs. ([7l)- ([TT]) . 
changes in the dielectric constant e of the microcavity 
result in the variations of the microcavity length Lc, for 
which the cavity photons and excitons in graphene are in 
the resonance, as well as in the changes in the polariton 
effective mass m and the polariton interaction strength g. 
It was found that the BEC density distributions for CuBr 
{e = 5.7), ZnSe {s = 9.1), CdTe {e = 10.2) and ZnO 
{e = 10.8) microcavities are qualitatively similar to those 
obtained for the microcavity synthesized from GaAs and 
considered above in this Section. To characterize the de- 
tails of the changes in the polariton BEC distribution in 
those microcavities, we plot in Fig. [6] the dependence of 
the dimensionless density of the condensate, n(r)Ax^, 



at the point ri = (2000 Ax, 1500 Ax), similarly to that 
shown in Fig. [5] for a GaAs-based microcavity. It is seen 
in Fig. [6] that the polariton density at the distance ~ 75 
/im from the excitation spot center gradually decreases 
with the decrease of the dielectric constant of the mi- 
crocavity material. However, this variation is slow: the 
change of the dielectric constant from s = 13 (GaAs) to 
5.7 (CuBr) only results in ~ 19% decrease of the BEC 
density from n ^ O.ISSA^-^ to n ^ O.lOSAx-^ Thus, 
the BEC density in a polariton trace in a microcavity 
with embedded graphene at large distances from the ex- 
citation spot is higher for the microcavity with high di- 
electric constant. 

In Ref. [S^ it was shown that the Kosterlitz-Thouless 
transition temperature Tc for polaritons in a microcavity 
with embedded graphene increases with the decrease of 
the microcavity dielectric constant e. On the other hand, 
from the consideration above it follows that the increase 
of e leads to higher BEC densities at a large distance 
from the excitation spot. These conclusions are consis- 
tent with each other because Ref. [3J describes the redis- 
tribution of the polaritons in the superfluid between the 
normal and superfluid components at finite temperature, 
whereas in the present work we consider the dynamics of 
the system at a given total number of the polaritons in 
the BEC in the low temperature limit T <C Tc. 

C. Effect of the polariton lifetime on the BEC 
dynamics 

Finally, to demonstrate the effect of the polariton life- 
time on the BEC dynamics, we studied the BEC spread- 
ing in a semiconductor quantum well in a microcavity for 
traditional, short-lifetime polaritons, for which we took 
r = 3 ps [9|. In Fig. [7] we compare the BEC density for 
the short-lifetime polaritons in a microcavity with GaAs 
QW with that for the long-lifetime polaritons in a micro- 
cavity with GaAs QW and with embedded graphene at 
^ = 0.5 eV. It is clearly seen that for the short-lifetime 
polaritons, even in the presence of an external force, the 
BEC is mostly located in the region where it is directly 
excited by the external pumping and its density rapidly 
decreases at large distances from the excitation spot. Fig. 
[3 also shows that the long- lifetime polariton BEC density 
in graphene with S = 0.5 eV is higher than that in the 
GaAs quantum weh in a microcavity at large distances 
from the excitation spot. 

D. Comparison of the BEC propagation in 

microcavity with an embedded quantum well and 

graphene 



From the results obtained in Sec. lIV Al - IIVCl it follows 
that the propagation dynamics of the long-lifetime polari- 
ton BEC in a microcavity with an embedded GaAs quan- 
tum well and in a microcavity formed by GaAs, ZnSe, 
CdTe, and CuBr with an embedded graphene layer are 
quantitatively similar to each other. However, the po- 
lariton BEC density in the trace for a GaAs-based micro- 
cavity with embedded graphene at (5 = 0.5 eV is higher 
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FIG. 7: Density of a long-lifetime polariton BEG in graphene 
for Uo = 5 meV and ^ = 0.5 eV, in a GaAs-based microcavity 
for t/o = 3 meV, and of short- lifetime polariton BEG in a 
GaAs-based microcavity with r = 3 ps for [/q = 5 meV. The 
force is Fo = 13 meV/mm. 

than that for a microcavity with GaAs QW, as is demon- 
strated in Fig. [71 The EEC density in the microcavity 
with embedded graphene sharply increases with the rise 
of the gap energy (5, as shown in Fig. [SJd. This allows 
one to utilize the gap energy in graphene as a param- 
eter that controls the polariton EEC propagation in a 
microcavity. The polariton EEC density far from the ex- 
citation spot gradually decreases with the decrease of the 
dielectric constant of the microcavity material, thus semi- 
conductors with higher dielectric constant provide better 
conditions for the observation of the polariton EEC prop- 
agation. 

V. CONCLUSION 

Through simulations of the exciton polariton EEC dy- 
namics by using the non equilibrium Gross-Pitaevskii 



equation, we demonstrated that the long-lifetime (r ~ 
100 ps) polaritons in a wedge-shaped microcavity can 
propagate over a macroscopically long distance ^ 500 
/im. This distance is large compared to that for short- 
lifetime polaritons in traditional experiments ^, i9|, ll6|, ll9|, 
|48J . The maximum of the polariton density in the EEC 
is shifted from the center of the excitation spot in the di- 
rection of the external force due to the exciton-polariton 
interaction. 

We also proposed to observe a polariton EEC propaga- 
tion in a microcavity with an embedded gapped graphene 
layer. It was found that the EEC density at large dis- 
tances from the excitation spot in a semiconductor quan- 
tum well and in gapped graphene are comparable with 
each other. However, in graphene there is an additional 
parameter that controls the long-lifetime polariton prop- 
agation, which is the energy of a gap in the electron and 
hole energetic spectra. The obtained results can be use- 
ful for practical applications of coherent polariton flow in 
a high-quality microcavity in working elements of inte- 
grated optical circuits 0-[S| • The advantage of graphene 
in a microcavity is that the propagation of a polariton 
EEC is dynamically tunable via electrostatic gating. 
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